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We study the ground state properties of doubly closed shell nuclei 16 O and 
40 Ca in the framework of Correlated Basis Function theory using state depen- 
dent correlations, with central and tensor components. The realistic Argonne 



^ ■ U14 and v' g two-nucleon potentials and three-nucleon potentials of the Urbana 



class have been adopted. By means of the Fermi Hypernetted Chain integral 
equations, in conjunction with the Single Operator Chain approximation, we 
evaluate the ground state energy, one- and two-body densities and electro- 
magnetic and spin static responses for both nuclei. In 16 we compare our 
results with the available Monte Carlo and Coupled Cluster ones and find a 
satisfying agreement. As in the nuclear matter case with similar interactions 

and wave functions, the nuclei result under-bound by 2-3 MeV/A. 
21.60. Gx, 21.10.Dr, 27.20. +n, 27.40. +z 



Typeset using REVTpjX 



I. INTRODUCTION 



The attempt to describe all nuclei starting from the same nucleon-nucleon interaction 
which reproduces the properties of two-, and possibly three-, nucleon systems is slowly 
obtaining its first successes. A set of techniques to exactly solve the Schrodinger equation 
in the 3<A<8 nuclei is now available: Faddeev 0, Correlated Hyperspherical Harmonics 
Expansion [0, Quantum Monte Carlo 0. Their straightforward extension to medium- heavy 
nuclei is however not yet feasible, both for computational and theoretical reasons. 

The Correlated Basis Function (CBF) theory is one of the most promising many-body 
tools currently under development to attack the problem of dealing with the complicate 
structure (short range repulsion and strong state dependence) of the nuclear interaction. 
The CBF has a long record of applications in condensed matter physics, as well as in liquid 
helium and electron systems. In nuclear physics the most extensive use of CBF has been 
done in infinite nuclear and neutron matter. The neutron stars structure described via 
the CBF based neutron matter equation of state is in nice agreement with the current 
observational data JU|. m nuclear matter CBF has been used not only to study ground 
state properties @,||J3] Du t also dynamical quantities, as electromagnetic responses ]8]|| and 
one-body Green's functions ||10|| . 

The CBF theory is based upon the variational principle, i.e. one searches for the mini- 
mum of the energy functional 

sm = [ I 1 > (i) 

in the Hilbert subspace of the correlated many-body wave functions 

V{1,2...A) = G(1,2...A)$(1,2...A), (2) 

where G(l, 2... A) is a many-body correlation operator acting on the mean field wave function 
$(1,2. ..A) (we will take a Slater determinant of single particle wave functions, (f) a (i)). In 
realistic nuclear matter calculations, the correlation operator is given by a symmetrized 
product of two-body correlation operators, Fij, 



(3) 

In principle richer choices for the operator (|3|) can be made by including explicit three-, 
or more- , nucleon correlations, which cannot be described by the product of two-body 
correlations. It is essential, however, that the two-body correlation Fij has an operatorial 
dependence analogous to that of the modern nucleon-nucleon interactions |TT|-[l5|. Nowa- 
days, sophisticated CBF calculations consider of the form: 

Fii = E (4) 

p=l,8 

where the involved operators are: 

Off 1 ' 8 = [1, <T, ■ (Tj, Sij, (L • S)ij] <g> [1, Ti ■ Tj] (5) 

and Sij = (3 r»j • (Ti r»j • (Tj — (Ti ■ (Tj) is the tensor operator. The correlation functions f p (r), 
as well as the set of single particle wave functions, are fixed by the energy minimization 
procedure. 

A key point in applying CBF is the evaluation of the many-variables integrals necessary 
to calculate the energy functional (|I|). A direct approach consists in using Monte Carlo 
sampling techniques (Variational MonteCarlo, VMC) [|l^]. However, the required numerical 
effort is such that, for realistic interactions and correlations, VMC can be efficiently used 
only in light nuclei. Actually, a realistic calculation of the ground state of 16 has been 
done in Ref. [|17| by using the so called Cluster MonteCarlo (CMC) method. In CMC the 
terms related to the scalar part of the correlation (p = 1) are completely summed by VMC, 
whereas the remaining operatorial (p > 1) contributions are approximated by considering 
up to four- or five-body cluster terms. 

An alternative to the MonteCarlo methodology is provided by cluster expansions and the 
integral summation technique known as Fermi HyperNetted Chain (FHNC) ||18|| , particularly 
suited to treat heavy systems. By means of the FHNC equations it is possible to sum infinite 
classes of Mayer-like diagrams resulting from the cluster expansion of the expectation value 
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G(1,2...A) =S 



i<j 



of the hamiltonian, or of any other operator. FHNC has been widely applied to both finite 
and infinite systems with purely scalar (state independent, or Jastrow) correlations. 

The case of the state dependent Fij, needed in nuclear systems, is more troublesome 
since the non commutativity of the correlation operators prevents from the development 
of a complete FHNC theory for the correlated wave function of Eq.(^j). For this reason an 
approximated treatment of the operatorial correlations, called Single Operator Chain (SOC), 
has been developed The SOC approximation, together with a full FHNC treatment 

of the Jastrow part of the correlation, provides an accurate description of infinite nucleonic 
matter ||. It is therefore believed that FHNC/SOC effectively includes the contribution of 
many-body correlated clusters at all orders. The evaluation of additional classes of diagrams 
in nuclear matter has set the estimated accuracy for the ground state energy to less than 1 
MeV at saturation density (p nm =0.16 fm~ 3 ) [p],pO[ . 



In a series of papers [f2T] -|23|| we have extended the FHNC scheme to describe the ground 
state of doubly closed shell nuclei, from 4 He to 208 Pb, with semi-realistic, central interactions 
and two-body correlations, either of the Jastrow type or depending, at most, on the third 
components of the isospins of the correlated nucleons. In Ref. we used FHNC/SOC to 
evaluate energies and densities of the 16 and 40 Ca nuclei, having doubly closed shells in 
the Is coupling scheme, with potentials and correlations containing operator terms up to the 
tensor components. In the 16 nucleus the comparison of our results with those of a CMC 
calculation confirmed the accuracy of the FHNC/SOC approximation estimated in nuclear 
matter. 

The present work is the extension of that of Ref. p4j| . The ground state properties 



of the 16 and 40 Ca nuclei are calculated within the FHNC/SOC formalism by using a 
complete, realistic nucleon-nucleon potential, with p > 6 components, and by considering 
also three-nucleon interactions. The two-nucleon interactions we have employed are the 
Argonne Vu |25| potential and the v' s reduction of the Argonne v±s potential. For 
the three-nucleon interaction we have adopted the Urbana models, Urbana VII [^(J with 
Argonne vu and Urbana IX || with Argonne v' 8 . In addition to the energy and the densities, 



we have also evaluated the static responses. They are the non energy weighted sums of the 
inclusive dynamical responses of the nucleus to external probes. We have studied the density, 
the electromagnetic and the spin static responses, both in the longitudinal and transverse 
channels. 

The paper is organized as follows: in section 2 we briefly present the interaction and 
the correlated wave function properties and recall the basic features of FHNC/SOC; section 
3 deals in short with the insertion of the spin-orbit components and of the three-nucleon 
potential; in section 4 we show and discuss the results for the energy, one- and two-body 
densities and static responses; the conclusions are drawn in section 5. 



II. INTERACTION, CORRELATED WAVE FUNCTION AND CLUSTER 

EXPANSION 

We work in the framework of the non relativistic description of the nucleus and use a 
hamiltonian of the form: 

# = ^Ev? + E%+ E %*■ (6) 

Zm i i<j i<j<k 

Very high quality phase-shift analyses of a large body of pp and np data have been recently 
carried out 0,0] • Building on this accurate data base, several nucleon-nucleon (NN) 



potentials have been constructed, like the updated Nijmegen interaction ||13|| , the CD Bonn 
interaction [pj|] and the Argonne t>ig (A18) 0] interaction. All of them include charge 
symmetry breaking terms in order to provide a precise fit to both the pp and np data. 

The structure of v% and at large inter-particle distances is dictated by meson ex- 
change processes. The long-range part of the NN interactions is determined by the one-pion 
exchange (OPE): 



= [YMi)°i ■ °i + T^Sij] r t ■ Tj . (7) 



where m n is the pion mass (138.03 MeV), (/^attv/^ 71 ") — 0.081, and Y n (r) and T v (r) are the 
Yukawa and tensor Yukawa functions, 



p— A" 

fir 



(8) 



-ftr 



//?' 



i 3 3 
1 + — + 



(9) 



/ir (/w) 2 

with /i ~ 0.7 fm -1 . The Argonne potentials simulate the effects of p exchange, not explicitly 
included | 25| , by modifying the Yukawa functions with a short-range cutoff, Y^ir) —> Y(r) = 
Y^(r)F cut (r) and T w (r) -> T(r) = T^(r)F 2 ut (r) withF cni (r) = l-exp(-cr 2 ) (c=2f m - 2 ). The 
intermediate and short-range parts of this class of potentials are mostly phenomenological 
and A18 is parametrized according to the operatorial structure: 



v 



j- E '■'•{n j )o>; r 

p=l,18 



where the first 14 components, 



o p ~- 



p=l,14 



1, 0"i ■ (Tj, Sij, (L ■ S)ij, L , L CTj • (Tj, (L • S)j • (g) [1, Tj • Tj] , 



give the isoscalar part, defining a t>i 4 -like potential (A14), 



(10) 



(11) 



v u,ij = E vP ( r ij)°ij- 

p=l,14 



(12) 



The remaining 4 components of A18 are of the isovector (r ijZ + Tj yZ ) and isotensor (Sri jZ Tj tZ — 
T i • Tj) type. The Argonne v' 8 potential (A8') is an eight operators reduction of A18 built 
to reproduce the isoscalar part of the full interaction in the S, P and 3 Di waves and the 
3 D 1 - 3 Si coupling 0. 



Other potentials use different parameterizations. For instance, the Nijmegen model fl2 
employs p 2 operators instead of L 2 . 

The longest range part of the three-nucleon interactions (TNI) involves a two pions 



exchange with the intermediate excitation of a A |27 



,.2tt 
ijk 



eye 



A2w E ( {^ij> ^ifc}{ T i ' T j; T i ' T fc} + T [Xij, X ik ] [t \ ■ T j , T 'j • T fc] ) , 



(13) 



where 



X i:j = Yfa^tTi ■ (t :i + Tfc^Sij , (14) 

and the symbols [, ] and {, } indicate the commutator and anti-commutator operators re- 
spectively. 

The Urbana class of TNI introduces an additional, repulsive, spin and isospin inde- 
pendent, short range term, 

v* k = U Y.T\r l0 )T\r lk \ (15) 

eye 

which simulates dispersive effects when integrating out A degrees of freedom. The total 
Urbana TNI are then given by the sum of the two terms defined above, i>y* = vfj k + v^ k , 
and the A 27T and U Q parameters are adjusted to provide a good fit to the binding energies 
of few-body nuclei and nuclear matter. In the Al4+Urbana VII (A14+UVII) model the 
values A2n= -0.0333 MeV and Uq— 0.0038 MeV have been fixed by variational calculations. 
In the most recent Al8+Urbana IX (A18+UIX) model the values A 2n = -0.0293 MeV and 
Uq— 0.0048 MeV have been obtained with a Quantum MonteCarlo calculation for 3 H and 
a variational calculation for nuclear matter. The values of A 2n are in good agreement with 
those predicted by the pure two-pions exchange model (A 27T ~ -0.03 MeV). 

The truncated A8' NN potential was introduced in Ref. because its simpler parame- 
terization allowed a simplification of the numerically involved Quantum MonteCarlo calcu- 
lations. The contribution of the missing channels was perturbatively evaluated. One should 
remark, however, that A8' was found to give a slight over bind. For this reason, the strength 
of the repulsive part of the TNI Urbana IX model was increased by ~30% to reproduce the 
experimental energies. The results presented in this paper have been obtained with the 
A14+UVII and A8'+UIX models, where the Urbana IX interaction has been redefined as 
above. 

The many-body wave function (|2[) contains two ingredients: the correlation functions 
f p (r) and the single particle states forming the Slater determinant $(1, 2. ..A). In our calcu- 
lations we have used a f^-type correlation, so neglecting the spin-orbit components. This 
choice will be discussed and commented later in the paper. 
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The best variational choice of would be given by the free minimization of the 
FHNC/SOC energy and the solution of the corresponding Euler equations, 5E/5Fij = 0. 
This approach is not practicable in finite nuclear systems, so we use an effective correlation 
obtained by the minimization of the energy evaluated at the lowest order of the cluster ex- 
pansion, Elo- The two-body Euler equations are then solved under the healing conditions 
/ 1 (r > di) = 1, f p>l (r > dp) = and requiring that the first derivatives of these functions 
at r = dp vanish. The healing distances d p are taken as variational parameters. In analogy 
with the nuclear matter case, in our calculations we adopt only two healing distances: d c for 
the four central channels and dt for the two tensor ones. Additional variational parameters 
are the quenching factors, a p , of the NN potentials in the Euler equations. More details are 
given in Ref. [|TJ| for nuclear matter and in Ref. |24| for finite nuclei. In the CMC calcula- 



tion of Ref. |T7j a nuclear matter Euler correlation was used and the nuclear matter Fermi 
momentum, kp , was used as a variational parameter. We follow here the same strategy. 

The single particle wave functions, 4> a (i), used in this work have been obtained either by 
solving the single particle Schrodinger equation with Woods-Saxon potential, 

Vws(r) = — W V ° w v (16) 

1 + exp[(r - Ro)/a \ 



or with a harmonic oscillator well, Vho{ t )i with oscillator length buo — yh/moo. The 
parameters V , _R , and a of Vws( r ) an d buo of Vno( r ) are also variationally determined. 

It is possible to express the expectation values of n-body operators in terms of n-body 
density matrices, and related quantities. In particular, the one- and two-body densities, 
P\{v\) and p2( r i) r 2) , defined as: 

Pi(r) = (E<Kr-r 4 )> (17) 

i 

and 

pf(r,r') = (E^-r J )5(r'-r J )Of J ), (18) 

are needed to compute the energy mean value (|I|). In the above expressions we indicate 
the mean value of an operator Q as: (Q) = (^\Q\^)/(^\^). Cluster expansion and Fermi 
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HyperNetted Chain theory provide a viable way to evaluate the densities both in infinite and 
finite Fermi systems. We present here only some of the basic features of the FHNC/SOC 
computational scheme. More complete discussions of the FHNC theory are found in Refs. 
fl8| , |2T|j for scalar correlations and in Refs. fli9| , |24|| for the operatorial case. 



In FHNC theory, with scalar correlations, the densities are written in terms of the cor- 
relation, /^(r), and of the nodal and elementary functions, iV^(ri,r 2 ) and El y (r\, r 2 ), rep- 
resenting the sums of the corresponding nodal and elementary cluster diagrams. The x(y) 
label characterizes the exchange pattern at the external points 1(2) and can indicate: a 
direct link, (d), if the particle is not exchanged and the point is reached by a dynamical 
correlation, h(r) = [/ 1 (r)] 2 — 1; an exchange link, (e), if the particle belongs to a closed 
exchange loop; a cyclic link, (c), if the particle belongs to an open exchange loop. The pos- 
sible combinations are: dd, de, ed, ee and cc. In infinite, homogeneous systems, as nuclear 
matter, the FHNC functions depend only on the inter-particle distance, ri 2 . The FHNC 
integral equations allow for computing the nodal functions, once the elementary ones are 
known. 

With the introduction of operatorial correlations, the nodal functions gain a state depen- 
dence, iVI^ 1 (i"i, r 2 ). The correlation operators do not commute anymore among themselves, 
[Fu, Fi 3 ] 7^ 0, and one must take into account the various possible orderings in fl3|). For this 
reason a complete FHNC treatment for the full, state dependent densities is, at present, not 
possible. The Single Operator Chain approximation was introduced in Ref. [19) for nuclear 
matter, and extended to Is doubly closed shell nuclei in Ref. [p4| , for f 6 correlations. The 



SOC scheme consists in summing those p > 1 chains, where each link may contain only one 
operatorial element and scalar dressings at all orders. Notice that operatorial dependence 
comes also on account of the exchanges of two or more nucleons, since the space exchange 
operator is given by P i<7 - = — Yl, v =c,a,T,<jT^ijl^- Here, and in the following, we may refer to 
the operatorial channels as c (p — 1), a (spin), t (tensor) and b (spin-orbit). The isospin 
channels have an extra r-label. 

The elementary functions, E xy (ri, r 2 ), represent an input to the FHNC equations, as 



they cannot be calculated in a closed form. The FHNC/0 approximation consists in ne- 
glecting all the elementary contributions. This seemingly crude approximation is actually 
based on the fact that the elementary diagrams are highly connected and have, at least, a 
quadratic dependence on the density of the system. These diagrams are not expected to 
produce relevant contributions in the relatively low density nuclear systems, whereas they 
are important in denser systems, like atomic liquid Helium. A test of the validity of the 
FHNC/0 approximation, and, in general, of the importance of the elementary diagrams, is 



provided by the degree of accuracy in satisfying the densities normalizations. In Ref. f24 
particular attention has been paid to the normalization of the one-body density, 



A = J rf 3 r lPl (n) (19) 



and to that of the central and isospin two-body densities, 

1 



1 = TFT in / d ^ / <*%^(ri, r 2 ), (20) 



-l = ^/^i |rf 3 r 2 p5( ri ,r 2 ). (21) 

Deviations of the sum rules from their exact values are due to {%) the approximate evalu- 
ation of the elementary diagrams and (ii) the SOC approximation. The first item has been 



already investigated in Ref. ||21|| , where it has been found that the most relevant corrections 
to the FHNC/0 sum rules come from the E^ ch diagrams, i.e. ee-elementary diagrams whose 
external points belong to the same exchange loop. These diagrams mainly contribute to the 
isospin saturation sum rule ( pT| ) and to the potential energy, if the interaction has large ex- 
change terms. This is understood if we notice that a four-point elementary diagram, linear 
in the central link h(r), belongs to E^ ch , as well as diagrams linear in the operatorial link, 



In Ref. |23| it has been shown that the one-body density sum rule is violated in 
FHNC/SOC by less than 1% and the two-body density normalizations ( pO| , pTD by ~9% 
in the worst case ( 40 Ca with tensor correlations). This is the same degree of accuracy found 
in Ref. @] for nuclear matter calculations. 

10 



In Ref. pi] the comparison of the 16 FHNC energies calculated for a purely Jastrow 



correlation with the exact VMC estimates shows an excellent agreement, with a difference 
of ~1% in the kinetic energy ( 24.61 MeV in FHNC vs. 24.33±0.21 in VMC) and of ~2% 
in the potential energy ( -22.07 MeV in FHNC vs. -21.56±0.25 in VMC). The 16 FHNC 
energies calculated with a f§ correlation and using the Vq truncation of the Urbana Vu 
potential |29| , have been compared with the results of a fifth order CMC calculation, which 



appeared to have reached a satisfying convergence. In this case we found a difference of less 
than 5% for the kinetic energy ( 31.16 MeV in FHNC/SOC vs. 29.45±0.33 in CMC) and 
of -7% for the potential one ( -35.47 MeV in FHNC/SOC vs. -33.03±0.31 in CMC). The 
FHNC/SOC calculation gives a binding energy per nucleon of -4.33 MeV to be compared 
with the -4.59±0.10 MeV value obtained by CMC. This difference is compatible with the 
nuclear matter estimates. 

III. SPIN ORBIT AND THREE BODY FORCES 

As already mentioned in the introduction, the novelty of this work with respect to Ref. 
2"4| is the inclusion of the spin-orbit and three-body terms of the potential. In this section 



we briefly show how to extend the FHNC/SOC formalism to consider these parts of the 
interaction. 

A. Spin Orbit potential 

The treatment of the spin-orbit interaction within the FHNC/SOC formalism has been 



discussed in details in Ref. |$0| for the nuclear matter case. In that paper, the evaluation of 
the mean value of the spin-orbit terms of the interaction , (v 7=b ' 8=bT ) , was done using a / 8 
correlation factor, so including all the 8 operatorial components in Eq.(^). 

Here we extend the nuclear matter formalism to the finite nuclei case for a f§ correlation. 
A correlation containing spin-orbit components would, probably, be more efficient from the 
variational point of view. However, these correlation terms introduce extra uncertainties in 
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the cluster expansion, and we have chosen to work with a simpler correlation in favor of a 
safer convergence. We shall discuss further this point and give some estimates of the f p>6 
corrections. 

In the FHNC/SOC scheme, the ground state matrix element, W, of a generic two-body 
operator, W, is split into four parts 

W = W + W S + W C + W CS . (22) 

where W$ indicates the sum of the diagrams with central chains between the fully correlated 
interacting points, connected by W, W s the sum of the diagrams having operatorial vertex 
corrections (or Single Operator Rings) touching the interacting points and central chains, 
W c the diagrams with one operatoral chain, SOC, between the interacting points, and, 
finally, W cs the diagrams containing both operatorial vertex corrections and chains. A more 
complete discussion on this topic is found in Refs. |T9|j2^] . 

The nuclear matter calculations of Ref. |30| show that two-body clusters provide the 
leading contribution to (v b ' br ). Three body separable diagrams, contributing to W s , gave 
the main many-body contributions and were of the order of 10% of the two body ones. 
Chain contributions were even smaller. These results have been confirmed in Ref. , where 
the authors found (v b+bT ) two-body / 'A =-2.29 MeV and a many-body contribution (W s + W c ) 
of -0.28 MeV at the nuclear matter saturation density, with the A18 potential. 

Relying on these facts, for the spin-orbit terms of the interaction, we have calculated 
only the Wq contribution of Eq.Q2"2]) with a fa correlation function. It turns out that, in this 
case, only the tensor correlations contribute, and the full expression reads as: 



(v^) Wo = -9 J d 3 n J d 3 r 2 /^(r 12 )^^(r 12 )/^(r 12 )/i c (r 1; r 2 ) (23) 
{p c (r 1 )p c (r 2 )K i ^A k ^ - 8(A^( ri ,r 2 ) - p ( ri , r 2 )) 2 ^'^^^A^A^Q(r 1 )Q(r 2 )} . 

In the above equation a sum over repeated indexes is understood, (i T , j T , k T ,l T ,m T = c, r), 
the coefficients K^ k , A 1 and A* are given in Ref. (TP| and the other FHNC functions are 



defined in 24 . 



12 



For the 16 O nucleus, we could check the reliability of this approximation against some 



CMC results ||31|| . Using a nuclear matter correlation and the A14 potential, we find 
(v b+bT )w /A = 0.56 MeV whereas CMC, with the same correlation, gives an extrapolated 
value of 0.62 MeV. We did not find the same kind of agreement for a fs correlation, where 
the Wq result is -1.23 MeV and CMC gives 0.03 MeV. This supports our choice in favor of 
the use of a simpler f$ correlation. 



B. Three body potential 

The evaluation of the mean value of the three-body interaction, = {vfj k ) + (v^ k ) 

closely follows the method developed in Ref. for nuclear matter. 

Diagrams 2.1-3 of the reference were shown to provide the relevant contributions to 
(vff k ), and diagrams 3.1 and 3.2 to (v^ k ). These diagrams, which we show in Figure [l], are 
those considered in our calculations. 

As an example, we show here how the nuclear matter expression of diagram 2.1 is ex- 
tended to our case. The explicit expressions for the remaining terms are given in the Ap- 
pendix. 

In diagram 2.1 the pairs of nucleon connected by the operators Xij (pairs 12 and 13 in 
the figure) are dressed at all orders by Jastrow correlations, whereas the remaining pair (23) 
bears the full operatorial correlations. Only the anticommutator part of vfj k contributes 
to {vfj k ) 2 .i- Spin-ispospin trace and spatial integration over nucleon 1 generate an effective 
two-body potential, acting on the 23 pair, having or and tr components only and depending 
on the exchange patterns, [xy] , of particles 2 and 3, 



with 



< /,fc (r 2 ,r 3 ) 



< / (r 2 ,r 3 ) 



E > 

k=(TT,tT 



,,eff,k, 
J xy ' 



r 2 ,r 3 )0 



k 

23) 



AA 2w J d 3 ri 



E 9 C XX >(*2, r i )V; c v (r 1 )^(r 1 , r 3 ) 



(24) 



(25) 
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where k = k a r, j a , k a = c, a, t), g c xy and V£ y are central partial distribution functions and 
vertex corrections defined in [pi], £ 2 { 3 are angular couplings given in Ref. |19] and X c (r) = 0, 



X a (r) = Y(r) and X*(r) = T(r). The allowed combinations for (x'y') are: dd, de, ed, and 
cc. In the last case, also (xy) should be (cc). 

The full expression for the mean value of the 2.1 diagram, (vfj c k ) 2 .i, is then given by: 

W*>2.1= \ J d 3 r 2 J d 3 r 3 P(r 23 )P(r 23 )h c (r 2 , r 3 )K^A^ (26) 
{v e // k (r 2 , r 3 ) K(r 2 )^(r 3 ) + p?(r 2 )iV d c e (r 2 , r 3 )C rf (r 3 ) + C d (r 2 )iV e c d (r 2 , r 3 )p<(r 3 ) 
+C d (r 2 ) (iV d c e (r 2 , r 3 )iV e c d (r 2 , r 3 ) + iV e c e (r 2 , r 3 )) C d (r 3 )] 
+^ e //,fc (r 2 , r 3 ) K(r 2 ) + C,(r 2 )iV e c d (r 2 , r 3 )] C d (r 3 ) 
+<P(r 2 , r 3 )C d (r 2 ) ^(r s ) + iV d c e (r 2 , r 3 )C d (r 3 )] 
+^P(r 2 ,r 3 )C d (r 2 )C d (r 3 )} 

-2 1 d\ 2 J d 3 r 3 f (r 23 )f (r 23 )h c (r 2 ,r 3 )C d (r 2 )C d (r 3 )K M K l ^ m A m A m 

{v e // k (r 2 , r 3 ) [iV c c c (r 2 , r 3 ) - p (r 2 , r 3 )] 2 + 2t#''*(r 2| r 3 ) [iV c c c (r 2 , r 3 ) - p (r 2 , r 3 )]} . 

Both commutator and anticommutator contributions are present in (f|^) 2 . 2 , while only 
the commutator part contributes to {vfj k ) 2 , 3 . 

In Table | we compare the TNI mean values for 16 calculated by FHNC/SOC with the 



MonteCarlo estimates 311 obtained with the same wave function. We have used the Urbana 



V three-nucleon interaction |28[ and the v§ truncation of the Urbana t>i4 NN interaction. 
The correlations have been obtained by the two-body Euler equation and we used harmonic 
oscillator single particle states with b = 1.54 fm. The 40 Ca results are also shown. 

For the Jastrow correlation the MonteCarlo results have been obtained with a VMC 
calculation, while for the fa correlation the calculation is of the CMC type. In 16 O the 
FHNC/SOC results are in satisfactory agreement with the MonteCarlo ones for both the 
classes of correlations. The 40 Ca results may be compared with those in nuclear matter, 
where, at saturation density, the Urbana V model gives (t>|^)=-2.32 MeV and (t>^ fc )=3.35 
MeV [[28 1 . 
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The change of sign of (vf^.) between the Jastrow and the operatorial correlation is due 
to the fact that the attractive contribution comes most from the tensor component of the 
effective potential, v e J^ tT {j2^ r%), which does not contribute in absence of tensor correlations. 



IV. RESULTS 

The results presented in this section have been obtained with v$ type NN interactions, 
either by truncating the A14 potential or by using the A8' potential of Ref. 0, together 
with / 6 correlations and Urbana three-nucleon potentials. 

We have estimated the f p>6 corrections, as well as the contributions from the p > 8 
components of the potential, by Local Density Approximation (LDA). In practice, if we 
define as AE nm (p) the sum of these corrections in nuclear matter at density p, we evaluate 
their contribution in the finite nucleus as: 

AE = \J rf3r Pi(r)A£ nm [ Pl (r)]. (27) 
We have already studied the accuracy of FHNC/SOC computational scheme against the 



results of CMC in Ref. [24]. In the present article we compare again the 16 results with 
those of CMC in order to test the accuracy in the calculation of the interaction terms we 
have added. For this reason we have computed the ground state energy for the A14+UVII 
model using a Jq correlation derived from the nuclear matter two-body Euler equations. 
The parameters of the correlation are: the Fermi momentum kp, and the healing distances 
d c , used for all the central channels, and dt for the tensor channels. Additional variational 
parameters are the quenching factors (3 P of the potential (see eq. 3.2 of |[9]]), and as in Ref. 
|2~4| we have taken f3\ = 1 and (3 P> \ = as- We have used the same set of single particle 
states, produced by a Woods-Saxon plus wine-bottle mean field potential, and correlation 
parameters of Ref. ||17|. However, small differences in the correlation come on account of the 



fact that our nuclear matter Euler equation does not contain the v p>6 components, contrary 
to Ref. IH7II. 
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The results of these calculations are shown in Table |H]. The Table gives also the cor- 
responding CMC results, as extracted from the reference, by subtracting, when reported, 
the three-body and spin-orbit correlations contributions. The AEqmc line gives the CMC 
corrections coming from (/, t>) p>6 and should be compared with the AE + (vj-s^j) sum (-0.5 
vs. -0.80 MeV/A) in the FHNC/SOC column. The ground state energy, E gs , is given by 
E gs = E v — T cm , where T cm is the center of mass kinetic energy. The Coulomb potential 
energy, (vcoui), has been included in E gs The root mean square radii, rms, are also reported. 

The energy found in CMC is -7.7 MeV/A, and it contains a -0.85 MeV/A contribution 
from explicit three-body correlations. Therefore, the FHNC/SOC result (-5.97 MeV/A) 
should be compared with an estimated CMC value of -6.85 MeV/A. If the CMC expansion 
has reached a satisfactory convergence, the discrepancies between the two calculations are 
due to the truncation in the FHNC/SOC scheme and to the aforementioned differences in 
the correlation. As it was found in Ref. [p4}| , the FHNC ground state energy differs from 



the CMC estimate by ~1 MeV/A, compatible with the estimated accuracy of the method 
in nuclear matter. In this calculation the AE corrections are small, somehow justifying a 
posteriori the use of LDA. It is remarkable that in the coupled cluster approach of Ref. 
the authors find E gs (CC)=-6.1 MeV/A with the two-body A14 model, close to our 



estimated £' 9S (FHNC/SOC)=-6.04 MeV/A (obtained without < v ijk > and < v Cou i >)• 

The nuclear matter energies per nucleon calculated with several interactions are given in 
Tab. |ffl]as a function of the density. The A18+UIX (A18) and A8'+UIX (A8') interactions 
provide energy minima at the empirical value of the saturation density. The results of 
the A18 column differ from those of Ref. because we have subtracted a perturbative 
correction, related to additional state dependence of the correlation. The second column 
of the Table shows the results for the A8'+UIX interaction with a correlation. The 
contribution of the f p>6 terms is given in the third column. For completeness, we also give 
the results obtained with the full A14 potential plus the UVII three-body force and the 
corresponding AE/A. The A14+UVII (A14) minimum of the equation of state is located at 
a density slightly higher than the empirical one. 
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The ground state energies of 16 and 40 Ca calculated with the A8'+UIX and the trun- 
cated A14+UVII interactions are shown in Table IV. We have adopted the nuclear matter 
correlation, where the nuclear matter density is used as a variational parameter (it means 
that, for a given p nm , we use the correlation function parameters minimizing the nuclear 
matter energy at that density). In addition, the energy has been minimized over the single 
particle potential (HO or WS) parameters. The Table also contains the kinetic energies 
computed with only the mean field wave functions. They are about half of the total kinetic 
energies obtained in the full calculations and this difference has to be ascribed to the corre- 
lations. This type of behavior is also found in nuclear matter. For example, at saturation 
density and with the A18+UIX hamiltonian, the nuclear matter calculations of Ref. || pro- 
vide a total kinetic energy of 42.27 MeV/A to be compared with the corresponding Fermi 
gas value of 22.11 MeV/A. As expected, the minimization with the WS potentials produces 
lower minima than those obtained with HO potentials. The comparison between the WS 
calculations done with the two different interactions shows small differences in the results, 
of the order of 7% in E gs . From the relative point of view the largest differences are those 
related to AE, which for A14 are 50% smaller than those of A8'. These rather similar values 
of the energies have been obtained with two quite different sets of single particle wave func- 
tions. This can be deduced by looking at the values of the WS parameters in Tab. [TV], and 
even better from Figure [| where the FHNC/SOC charge distributions are compared with 
the IPM ones and with the empirical densities taken from the compilation of Ref. |I3] . The 



charge densities have been obtained by folding the proton densities with the electromagnetic 
form factors of Ref. |34] . 



We do not obtain a satisfactory agreement between the computed densities and the 
empirical ones. However, we remark that, for a given type of single particle wavef unctions, 
either HO or WS, we find shallow energy minima with respect to variations of the mean field 
parameters around the minimum itself. This may indicate that charge distributions and rms 
radii are sensitive to details of the many body wave function which have small effects on the 
energy calculation. To better study this aspect, we have done calculations by using a set of 
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single particle wave functions fixed to reproduce at best the empirical charge densities. In this 
case, the only variational parameters are those related to the correlation functions. Actually, 
we have varied only d t and kept fixed its ratio with d c . The results of these calculations 
with the A8'+UIX interaction are shown in Tab. [V] and in Figure [3| We observe that a 
large change in the single particle wave functions produces small variations in the energy 
values, 1.3% in 16 and 4.8% in 40 Ca, within the accuracy of the FHNC/SOC scheme. The 
analysis of these results shows that the values of the kinetic and potential energies have 
considerably changed, while their sum remains almost constant. The agreement with the 
empirical densities has clearly improved, as one can also see in Figure f| where the elastic 



electron scattering cross sections calculated in Distorted Wave Born Approximation |35 
with the FHNC/SOC charge densities are compared with the experimental data [|3~6|1 . The 
best agreement with the data is produced by the densities of Figure [3[ However, also in this 
case the high momentum and large energy data are not well described in both nuclei. 

From the comparison between the dashed and the thin full lines in Figures ^| and |^ we can 
inspect the effects of the correlation on the charge densities. In 16 the correlations enhance 
the densities with respect to the IPM ones without substantially changing the shape of the 
oscillation. The effect of the correlations on the 40 Ca densities is negligible. These results 
confirm the findings of Ref. [37j where the charge densities have been calculated up to the 
first order in the correlation lines. On the other hand, in Ref. [17] the correlations produce 
much larger deviations from the IPM densities. The reason of the different effect of the 
correlations on the densities of 40 Ca and 16 is presently object of investigation. 

Short range correlations strongly affect the two-body densities fll8l). Their effect is 
evident in the two-nucleon distribution function, p 2 (ri 2 ), defined as: 

P2(ri 2 ) = jJ d 3 R 12 p c 2 {r 1 ,r 2 ) (28) 

where R42 = \ {v\ + is the center of mass coordinate. In an analogous way, we define 
the proton-proton distribution function, p pp (r 12 ), as: 

P P p(r 12 ) = zf d 3 Ri 2 p pp (r 1 ,r 2 ), (29) 



where the pp-two-body density is: 

Mri,r 2 )= (ZSirx-v^-vj) (1±^) (30) 

Figure ^| shows ^2(^12) an d p PP (^i2) with the wave functions of Tab. [V], compared with 
the IPM densities. The reduction of the correlated distribution functions at small r 12 -values 
is due to the repulsive core of the interaction. We have calculated the FHNC/SOC 16 
distribution functions also for the A14+UVII model and we found that they are rather similar 
to the A8'+UIX ones {p™ ax (AU + UVII) ~0.081, p% ax (A8' + UIX) ~0.089, p™ ax (AU + 
UVII) ~0.032 and p^ x (A8' + UIX) ~0.032), and in good agreement with those computed 
in Ref . . 



The interest in calculating the two-body densities is also related to the possibility of using 
their Fourier transforms to analyze several integrated nuclear responses. The responses of 
the nucleus to external probes, either of electromagnetic or hadronic type, can be related to 
the dynamical structure functions, Sx(q,u), given by 

S x (q, uo) = K*/|Ojf |*o)r%i - w - uj), (31) 

where Ox is the operator producing the fluctuations around the ground state \l/o- In the 
above equation the sum runs over the intermediate ^/-states with energy ui. The non 
energy weighted sums of Sx(q,u) give the static structure functions (or, simply, structure 
functions, SF), Sx(q), as: 

Sx(q) = J S x (q,uj)dw = (^ \O x O x \^ ). (32) 

In the case of density fluctuations, the operator is: 

Pq = exp(«q- ri) (33) 

i=l,A 

and the lower limit cu-integration in fl3"2] ) is taken in an appropriate way to eliminate the 
contribution of the elastic scattering. The density SF, S(q), is then 
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S(q) = l + -f rf 3 r 1 rf 3 r 2 exp(zq-r 12 ){p^(r 1 ,r 2 ) - pi(ri)pi(r 2 )} . (34) 

From the normalizations of the one- and two-body densities, one obtains S(q = 0) =0. 
The response to charge fluctuations is driven by the operator 

p c ,q = J2 ex p(*q • r i) ( + Tl ' z ) > ( 35 ) 



1=1,^4 x " ' 

which is responsible also of the electromagnetic longitudinal response, if the small (~ 2%) 
contributions of the neutron magnetic moment and of the meson-exchange currents are 
disregarded. The longitudinal response is measured in inelastic electron-nucleus experiments 
and its energy integral gives the longitudinal SF, S^q), or Coulomb sum. In the nuclei we 
are studying, the explicit expression of S^q) reads: 

S L {q) = 1 + T^/ ^ 3 n^ 3 r2exp(zq-r 12 )|p^(r 1 ,r 2 ) + ip5(r 1 ,r 2 )-p 1 (r 1 )p 1 (r 2 )|. (36) 

In Figure |6] we present the density (upper panel) and charge (lower panel) SF for 16 0, 
calculated with the wave functions of Tab. M. The FHNC/SOC results obtained from our 
calculations are shown by the white triangles. The q = value is very sensitive even to 
small violations of the normalizations of the densities. For instance, the two-body density 
normalization of Eq . (|20"D , with the Tab. |V| parameters, is violated in 16 by only 5.3% 
and that of the one-body density, Eq. (|i~9D , by 2.0%. These acceptable normalization errors 
produce the large values S(q = 0) = 1.61 and Si(q = 0)=0.86. 

The SF obtained after properly renormalizing the densities are given in Figure |6| as 
full curves. The other curves show the IPM SF (dashed lines) and those obtained with 
only Jastrow correlations. The renormalization is effective only at small g-values and for 
q >1 fm -1 the SF remain unchanged. Analogous calculations done in 40 Ca show a similar 
behavior. Hence, the conclusions about the large-g importance of the correlations in the SF 
and densities are not affected by the normalization problems. In agreement with the findings 
of Ref. pTfl , our results show that the correlations, both of the Jastrow and operatorial types, 



lower the SF at large q. 
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In Figure [7| we compare the Coulomb sum rules calculated for 16 O and 40 Ca with the 
A8'+UIX interaction with the experimental estimates done by analyzing the set of world 
data on inclusive quasi-elastic electron scattering p8[ experiments in 12 C, 40 Ca, and 56 Fe. 



The Figure also shows the nuclear matter Coulomb sum for the A14+UVII model from 
Ref. JJIJ. The finite nuclei results are in complete agreement with the latest analysis of the 
experimental data, where a detailed study of the electron scattering world experiments and 
a proper inclusion of the large energy tail in the dynamical response were carried out. The 
nuclear matter results fail in reproducing the data at the lowest q values where finite size 
effects can be still important. 

Besides the density and charge SF also the isovector spin longitudinal and transverse (ISL 
and 1ST) SF are of experimental interest since they can be extracted from polarized proton 
and neutron scattering cross section data. Experiments of this type on 12 C and 40 Ca nuclei 
40| did not confirm the prediction of Random Phase Approximation [41 and Distorted- 



Wave Impulse Approximation [^] calculations of a large enhancement, with respect to unity, 
of the ratio of the ISL to the 1ST response at small energies. 

The ratio has been calculated for nuclear matter within the CBF theory, using a Jq 
correlation together with the Urbana t>i4 +TNI potential The computed average en- 
hancement was ~ 20%, compatible with the data in heavy nuclei at energies below 100 MeV. 
However, the nuclear matter calculation did not take into account the strong distortion of 
the emitted nucleon wave function. Variational and cluster Monte Carlo wave functions 
were used in Ref. |44] to evaluate the integrated spin responses in light nuclei and 16 0. A 



maximum 25% enhancement of the ratio was found in 16 0. 
The fluctuation operators in the isovector spin responses are 



PaL.q = eX P(«l ■ r i)(°"i ' q)( T i ' T ) ( 37 ) 

l = l,A 



in the longitudinal case, and 



PaT,d = ex P(^q • r i)(° r i x q)( T i • T ) ( 38 ) 

i=l,A 
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in the transverse one. In the above equations we have indicated with T a unit vector in 
the isospin space. If the nucleus has zero isospin, then the response does not depend on the 



direction of T, except for the small Coulomb effects. Following the treatment of Ref. [g4 
we obtain for the ISL structure function, S a L(q), the expression: 

S*L(q) = = 1 + / d3r ^ K(ri, r 2 )j (gr 12 ) - pf(v u r 2 )j 2 (qr l2 )} , (39) 

and for the 1ST one, S a T{q)'- 

SMq) = = 1 + <LJ d3rid3r2 { P7{ru r ^ qr ^ + ^^(ri, r 2 )j 2 (gr 12 )} . (40) 

The ISL and 1ST structure functions for 16 and 40 Ca, calculated with the wave functions 
of Tab. [V|, are shown in Figure |] as a function of the momentum transfer. For the 16 nucleus 
we show also the Jastrow results, identical for both ISL and 1ST. Central correlations do 
not differentiate between the longitudinal and transverse responses because of the lack of 
tensor correlations. Our 16 results are very close to those of Ref. Q| obtained with the 
A14+UVII interaction. The lowest panel of Figure shows the S a L(q) / S a T(q) ratio. We 
observe that the maximum enhancement is ~25% in 16 and ~38% in 40 Ca, close to the 
experimental value. 

V. CONCLUSIONS 



This paper is the natural extension of the work of Ref. [p4fl . We have added to the 



FHNC/SOC computational scheme for doubly closed shell nuclei in Is coupling, the con- 
tribution of the spin-orbit and three-body interactions. Our calculations of 16 and 40 Ca 
nuclei have been done considering a f§ type of correlation and v$ potentials: the A8'+UIX 
H and a truncated version of the A14+UVII [f25| potentials. The contribution from the 
remaining momentum dependent terms of the correlation and of the potential are estimated 
by means of a local density interpolation of their nuclear matter values. 

The main results of the paper are given in Tabs. [TV] and |V| where the 16 and 40 Ca 
ground state energies per nucleon are presented. Their values are ~ 2-3 MeV/A above the 
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experimental ones, consistently with the CBF results in nuclear matter. Additional lowering 
of the energies may be obtained by (i) three-body correlations, (ii) perturbative corrections 
to the two-body correlations. It has been already mentioned that three-body correlations 
have been found to provide an extra 0.8 MeV/A binding in 16 for the A14+UVII model. 
Perturbative corrections have been taken into account in nuclear matter either by the method 
discussed in Ref. || or by the inclusion of the second order two-particle two-hole contribu- 
tion in correlated basis perturbation theory |45| , 46|| . Both approaches lower the energy by 



AE2 ~2. MeV/A. The nuclear matter case gives a strong indication that the inclusion of 
these corrections should be pursued and that their quantitative consistency in finite nuclear 
systems needs to be numerically checked. 

A complete minimization over all the parameters of the wave functions, both in the 
correlation and in the mean field, has led to a marked disagreement between the CBF and 
the empirical charge densities in the low distance region. However, calculations with different 
sets of single particle wave functions, reproducing the empirical densities in IPM, provided 
energies differing from the best minima only by a few percent. The CBF scheme does not 
appear to be very sensitive to the details of the mean field basis in a parameters region 
around the variational minimum. Nevertheless, these details become relevant for a correct 
description of the one- and two-body densities. The introduction of additional constraints 
during the minimization process may be necessary in order to avoid these ambiguities. 

Our results show that the short-range correlations produce small effects on the charge 
density distributions, especially in 40 Ca where the FHNC scheme is supposed to perform 



better. These findings are in agreement with those of Ref. |37j where the same kind of 
nuclear matter correlations have been used. The sensitivity of the charge distributions 
to the state dependent short-range correlations requires further investigations to be fully 



clarified. In effect, the VMC one body densities of ]T7[ have larger dependence on these 
components. 

In addition to the ground state energies, we have studied the static structure functions. 
The FHNC results for the Coulomb sum rule fully agree with the empirical values. The 
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ratio between the ISL and 1ST SF shows an enhancement between 25% in 16 and 38% in 
40 Ca, in agreement with those of Refs. P31B3| and just slightly higher than the experimental 



estimates. 

From this work we can conclude that realistic variational calculations for medium-heavy, 
doubly-closed shell nuclei in Is coupling scheme with modern, sophisticated potentials are 
not only feasible, but have also reached the same degree of accuracy as in nuclear matter. 

There are several natural extensions we envisage. For instance, the inclusion of three- 
body correlations, the study of N^Z closed shell nuclei and the development of the 
FHNC/SOC formalism for the jj coupling scheme. Work along these directions is in progress. 
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VII. APPENDIX 

In this Appendix we give the explicit expressions of the remaining diagrams contributing 
to the three-nucleon potential expectation value in FHNC/SOC. 

<^>2. 2 = 6A 27T J d 3 ri Jd 5 r 2 J ^3 ; 9^ X ' CT ( ri2 ) J 9l^ XmCT(ri3) (41) 
■ R i,oMm a + SAH^AH^] p c 3 (r u r 2 , r 3 ), 

<^> 23 = 12A 2n J d» ri / d\ 2 J rf3r 3^^^ iCT (n 2 )^^X^(r 13 ) (42) 

" i*l*n„ _ K iMn„ A n^ g*™*"* pfa^ r 2 , r 3 ), 

( v 5khi = \J d?n J Sr 2 J dVauf^r!, r 2 , r 3 ), (43) 
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«3fc)3.2= ^ / d 3 r a / d 3 r 2 I d 3 r 3 v* 3 



1 

2 



f(r 23 )f k (r 23 ) 
\Hr 23 W 

PW r i' r 2 , r 3 )^ fc + P3.ex C h( r i ; r 2, r 3 )A r 



(44) 



The L* jfc and i?^ fc ' matrices are given in Ref. |19] and Ref. |2"Bfl , respectively. The central 
three-body density, pg(ri, r 2 , r 3 ), is written in superposition approximation as: 

p 3 (ri,r 2 ,r 3 ) = £ #L'( r i, ^)V^ ly/ (r 2 )g c y , y (r 2 , r 3 )V y c z (r 3 )g c zz ,(r 3 , n) V^ x (ri) (45) 

£yza;'?/,z'=<i,e 

- 8^(ri, r 2 )V; c c (r 2 )^ c (r 2 , r 3 )^ c c c (r 3 )^ c c c (r 3 , n)V£(ri). 



As already stated, only the ee combination is not allowed at a given vertex. The exchange 
and direct three-body densities, P 3>ea ;c/i/dir( r i) r 2? r 3); are given by those parts of the full 
density where nucleons 2 and 3 belong or not to the same exchange loop. 
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TABLES 

TABLE I. TNI expectation values per nucleoli, in MeV, in 16 and 40 Ca for the Urbana 
V model. The /i( 6 ) column gives the energies for the Jastrow (operatorial) correlation. The 
superscript C(A) indicates the commutator (anticommutator) contribution. 



/i(FHNC) /i(VMC) /e(SOC) /e(CMC) 







16q 






/ 2tt,C\ 
( V ijk ) 


-0.17 


-0.16 


-0.90 


-0.86 


1 2vr,A\ 
( V ijk ) 


0.74 


0.70 


-0.39 


-0.44 




1.33 


1.28 


1.65 


1.57 






40 Ca 






( V ijk ) 
{ V ijk ) 
(vfik) 


-0.24 
1.80 
2.75 




-1.50 
-0.26 
3.20 
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TABLE II. FHNC/SOC and CMC energies in Mev/A for 16 with the A14+UVII interaction 
and the CMC single particle potential and nuclear matter correlation, rms in fm. See text. 





FHNC/SOC 


CMC 


(T) - T cm 


37.23 


32.0 


(V6,ij) 


-42.48 


-38.2 


<W7-8,ij> 


-0.85 






-2.83 






1.90 




(%'fc) 


-0.93 


-1.1 




1.00 


0.9 


AE 1 


0.05 








-0.5 


E gs 


-5.97 


-6.9 


rms 


2.44 


2.43 



TABLE III. Nuclear matter energies per nucleon for A18+UIX (A18), A8'+UIX (A8') and 
A14+UVII(A14). The AE/A columns list the spin-orbit correlation and the p >8 potential com- 
ponents corrections. Densities in fm -3 , energies in MeV. 



Pnm 


E/A(A18) 


E/A (A8') 


AE/A (A8') 


E/A(A14) 


AE/A (A14) 


0.04 


-4.11 


-5.50 


-0.30 


-5.25 


-0.80 


0.08 


-7.46 


-8.45 


-0.82 


-8.43 


-0.38 


0.12 


-9.42 


-10.31 


-1.48 


-10.63 


0.00 


0.16 


-10.05 


-10.87 


-2.16 


-11.99 


0.61 


0.20 


-8.74 


-10.06 


-3.02 


-12.37 


0.27 


0.24 


-5.66 


-7.75 


-3.96 


-11.74 


0.50 
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TABLE IV. Energies in MeV/A for 16 and 40 Ca obtained with A8'+UIX+/ 6jnm and with 
truncated Al4+UVII+/ 6 . nm . In the upper part of the table we show the values the minimization 
parameters, bnoi a o> Ro, are expressed in fm, p nrn in fm -3 all the other quantitites in MeV/A. The 
energy experimental values are -7.97 MeV/A for 16 and -8.55 MeV/A for 40 Ca. The experimental 
rms radii are 2.73 and 3.48 fm for 16 and 40 Ca, respectively. 



0(HOW 16 0(WS) A8 , 16 0(WS) A14 40 Ca(HO) A8 , 40 Ca(WS) A8 , 40 Ca(WS) A14 



bno 


2.00 






2.10 






V 




-42.00 


-36.00 




-50.00 


-41.50 


Rq 




3.60 


3.80 




5.30 


5.00 


a 




0.55 


0.55 




0.53 


0.55 


Pnm 


0.09 


0.09 


0.09 


0.09 


0.09 


0.09 


(T) 


22.57 


27.34 


25.91 


30.02 


30.58 


30.90 


{T)ipm 


11.64 


13.82 


12.20 


14.09 


14.15 


14.10 


( v S,ij) 


-27.49 


-32.48 


-31.05 


-38.03 


-38.68 


-39.08 


(vfj k ) 


-1.15 


-1.66 


-1.35 


-1.94 


-1.94 


-1.74 


(vfjk) 


1.26 


1.98 


1.56 


2.46 


2.37 


2.11 




0.48 


0.58 


0.52 


0.18 


0.20 


0.19 


(VCoul) 


0.78 


0.86 


0.84 


1.85 


1.85 


1.87 


AE 


-0.66 


-0.94 


-0.50 


-0.95 


-0.96 


-0.70 


Egs 


-5.18 


-5.48 


-5.11 


-6.77 


-6.97 


-6.50 


rms 


3.03 


2.83 


2.93 


3.65 


3.66 


3.66 
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TABLE V. Energies in MeV/A for 16 and 40 Ca obtained with A8'+UIX+/ 6jnrn and with WS 
mean field potentials fixed to reproduce the empirical charge densities. Rq and ao are expressed in 
fm, p nm in fm -3 , and the other quantities in MeV/A. 





16 o(wsw 


40 Ca(WSW 


V 


-53.00 


-50.00 


Ro 


3.45 


4.60 


ao 


0.7 


0.5 


Pnm 


0.09 


0.09 


(T) 


32.64 


38.15 


(T)ipm 


15.35 


16.82 




-37.79 


-46.34 


& 

\ I J hi i 


-2.36 


-2.98 




3.00 


3.94 


E 

J -'cin 


0.64 


0.23 


{vCoul) 


0.94 


2.10 


AE 


-1.21 


-1.28 


Egs 


-5.41 


-6.64 


rms 


2.67 


3.39 
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FIGURES 

FIG. 1. Cluster diagrams considered for the three-body force expectation value. The 2.1 — 3 
diagrams are related to the {vfj k ) part of the force and the 3.1 — 2 diagrams are related to {vfj k ). 
The points denote the particle coordinates. The dashed, wavy and double-wavy lines denote 



generalized scalar, operator and single-operator ring correlation bonds, respectively. See Ref. [28] 
for more details. 

FIG. 2. FHNC/SOC charge densities related to the results of Tab. f[V| (thin full lines) 
compared with the IPM densities (dashed lines) and with the empirical ones (thick full lines). 



FIG. 3. The same as Figure || for the calculation of Tab. |V| 

FIG. 4. Electron scattering elastic cross sections for 16 (left) and 40 Ca (right). The full 
lines have been produced by the FHNC/SOC densities of Figure ||. The other two lines by the 
WS FHNC/SOC densities of Figure 0. The dashed lines correspond to the A8' densities while the 
dotted lines to the A14 ones. 

FIG. 5. Two-nucleon (upper panels) and proton-proton (lower panels) distribution functions 
for the A8'+UIX interaction. The dashed lines give the IPM results, the solid lines the FHNC/SOC 
ones. 

FIG. 6. Density (upper panel) and longitudinal (lower panel) structure functions for 16 (left) 
for the A8'+UIX interaction. The IPM (dashed lines), FHNC/SOC (solid lines), non-normalized 
FHNC/SOC (triangles) and Jastrow (stars) results are given. 

FIG. 7. Coulomb sums for 16 0, Ca and nuclear matter (n.m.) compared with the experi- 
mental results from 12 C, 40 Ca and 56 Fe. 

FIG. 8. Isovector spin longitudinal (solid) and transverse (dashed) structure functions for the 
16 and 40 Ca nuclei. For 16 we show also the Jastrow results (dotted line). The lowest panel 
gives the ratios of the isovector spin longitudinal and transverse structure functions for both nuclei. 
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